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SPHERICAL MATRIX ENSEMBLES 


GENE S. KOPP AND STEVEN J. MILLER 


Abstract. The spherical orthogonal, unitary, and symplectic ensembles (SOE/SUE/SSE) Sp{N, r) 
consist of N X N real symmetric, complex hermitian, and quaternionic self-adjoint matrices of 
Erobenius norm r, made into a probability space with the uniform measure on the sphere. Eor each 
of these ensembles, we determine the joint eigenvalue distribution for each N, and we prove the 
empirical spectral measures rapidly converge to the semicircular distribution as N oo. In the 
unitary case (/3 = 2), we also find an explicit formula for the empirical spectral density for each N. 
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1. Introduction 

After writing this paper, it was brought to our attention that the 
our spherical ensembles have been studied previously by physicists 
Caer and Delannay, under the name of fixed trace ensembles HCDI, IDCI1 . 
Delannay and Caer prove results equivalent to our Theorem 1.4 (on 
moments of spherical ensembles) for 13 = 1,2 in Appendix B of [[DC]], 
and our Theorem 1.5 (on empirical spectral density) in Section 3.1 
of [IDCII . Our methods differ from those of Delannay and Caer: We use 
the method of moments, and they use the Laplace transform method. 
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Random matrix theory has grown enormously in results and scope from Wishart’s HWisH 1928 


paper. Through the work of Wigner [ |Wigl[ |Wig2[ |Wig3[ |Wig4[ |Wig5| , Dyson [ |Dyl[ |Dy2[ |, and 
others, it now successfully models a variety of systems from the energy levels of heavy nuclei and 
zeros of L-functions HConl IPorU to bus routes HBBDSi IKrSell ; see for example HAGZi IFori IMehi 
IMT-BH for detailed expositions of the general theory, and IlFMi |HayP for a review of some of the 
historical development and applications. 

Much of the subject revolves around determining properties of the spectra of self-adjoint ma¬ 
trices. In this paper we concentrate on the density of states, though there has been significant 
progress in recent years on the spacings between normalized eigenvalues (see for example HERSYi 
lESYl ITvn ITV2II 1. The limiting spectral measure has been computed for many ensembles of pat¬ 
terned or structured matrices, including band, circulant, Hankel and Toeplitz matrices, as well as 
adjacency matrices of d-regular graphs jBBDSi IBasBoli IBasBo2[ IBanBoi IBCGl IBHSll IBHS2[ 
IBM [BDJl IGKMNl iHMl ITMPl iKaR IKKMSXl ILWI iMcKl iMel ISchll . Many of these papers show 
that when a classical ensemble is modified by enforcing particular, usually linear, relations among 
the entries, new limiting behavior emerges. We study spherical ensembles introduced in HMT-BII 
(see Research Project 16.4.7). While in general it is difficult to obtain closed-form expressions for 
spectral statistics of structured ensembles, we are able to exploit the symmetry of the spherical en¬ 
sembles to give explicit formulas for the joint spectral distribution and empirical spectral density. 
We prove a semicircular law and give numerical evidence of GUE spacing statistics. 

As a set, the N x N spherical unitary ensemble of radius r consists of x Hermitian 
matrices 


( 


A = 


fll2 


flu 

- bi2s/^ 


fll 2 + ^12\/~1 
0-22 


OlN + \ 

02N + b2Nsf-^ 


\ OiM — 6iAr\/^ a2N — &2Ar 


OnN 


( 1 . 1 ) 




satisfying the relation 


X] 4 + 2 ^ (a-j + bij) 


= r 


( 1 . 2 ) 


i<j 


or more succinctly, Tr(A^) = r^. This set of matrices becomes a probability space with the uniform 
measure on the (compact) ellipsoid defined by Tr(A^) = r^. One interpretation of the name is that 
Tr(A^) = is a sphere in the Erobenius norm (see below). 

The spherical ensemble may be understood to arise naturally from the problem of normaliza¬ 
tion. Consider a Hermitian matrix A chosen from a Gaussian unitary ensemble, with volume form 
g- Tr(A 2 )/ 2 g^^ j-Q ^ constant). To study the limiting behavior as the dimension N ^ oo, it’s stan¬ 
dard to normalize the eigenvalues by a factor of 1/ \/N. With this normalization, the eigenvalue 
density converges to Wigner’s semicircle. 


fwig{x) 


if |a:| < 1, 
0, otherwise. 


(1.3) 
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One could instead normalize the matrices themselves before looking at the eigenvalues. We can 
normalize a Hermitian matrix A by sending it to where | is the Frobenius norm (see (ll.SI) ). 

This “pre-normalized GUE” is just the spherical unitary ensemble of radius ^/JV- 

We also consider spherical orthogonal and symplectic ensembles. In all cases, we prove that 
the eigenvalue density converges to Wigner’s semicircle as the dimension N ^ oo. In the fixed¬ 
dimensional case, however, the empirical and joint spectral distributions exhibit new behavior. We 
provide explicit formulas and investigate their analytic properties. 

1.1. Definitions. Let F be a finite-dimensional real, complex, or quatemionicQ vector space with 
a Hermitian inner product (—, —). For any A G End(l/), let A* E End(l/) be the adjoint, so 
(Av, w) = {v, A*w). The Frobenius inner product on End(l/) is defined by 

{A,B) := Ai{AB*), (1.4) 

and induces the Frobenius norm 

1^1 := s/{A,A) = ^Tt{AA*). (1.5) 

When A is written as a matrix (with respect to an orthonormal basis on V), the Frobenius norm of 
A equals the square root of the sum of the norm squares of the entries of A. If A is self-adjoint, 
then I A| = also equals the square root of the sum of the squares of the eigenvalues of A. 

For our purposes, an iV x iV matrix ensemble is simply a subset of where K = M, C or 

HI, along with probability measure. For a matrix ensemble 8, the underlying set is also denoted by 
8, and the associated volume form is denoted by dsA. The ensembles we examine are all ensem¬ 
bles of self-adjoint matrices and thus have real eigenvalues, so we assume that in the discussion 
that follows. 

If it exists, the empirical spectral density of the ensemble 8 is the unique (up to a measure zero 
set of values) real-valued function f{x,8) such that 




g{x)f{x,8)dx 


( 1 . 6 ) 


for all g E T£M), where Ai(A) A 2 (A) < • • • < Xn{88) are the eigenvalues of A. The 
moment of f{x,8) is 


mt{£) := = 4 f'a{A’‘)dcA, 


N 

and its corresponding characteristic function is 


N 




k=0 


and is related to the density by 


fix, 8) = 


271 


k\ 


e-^^^fiit,8)dt. 


(1.7) 


( 1 . 8 ) 


(1.9) 


Let Sn be the symmetric group on{l,2,...,iV}. The joint spectral distribution is the measure 
on the set M^/ Sn of N unordered points induced by pushing forward the measure on 8 along the 


*A quaternionic vector space is a free (H, H)-bimodule, where 11 denotes the Hamilton quaternions. A finite¬ 
dimensional quaternionic vector space is necessarily isomorphic to some H". 








map A !-)■ (Ai(y4), A 2 (A),..., Ajv(A)). We consider it as a measure on by symmetrizing it and 
normalizing it to be a probability measure. 

We can now carefully state the ensembles of interest in this paper. 

Definition 1.1 (Gaussian ensembles). The Gaussian orthogonal / unitary / symplectic ensembles 
(GOE/GUE/GSE) are denoted q) and consist of N x N real symmetric, complex Hermitian, 

and quaternionic self-adjoint matrices, made into a probability space with the density 

dG,iN,q)A := Cp{N,q)e-f^^\^\"/^dA. ( 1 . 10 ) 

Here f3 = 1,2,4/or the GOE, GUE, and GSE, respectively, q is a positive real number, and 
q) is chosen to give a probability density. 

Definition 1.2 (Spherical ensembles). The spherical orthogonal / unitary / symplectic ensembles 
(SOE/SUE/SSE) Si 3 {N, r) consist of N x N real symmetric, complex hermitian, and quaternionic 
self-adjoint matrices of Frobenius norm r, made into a probability space with the uniform measure 
on the sphere. 

1.2. Results. For all spherical ensembles, we describe the joint eigenvalue distribution. 

Theorem 1.3. Let (3 G {1,2,4}. The joint eigenvalue distribution for Sp{N,r) (on vectors of 
eigenvalues A G M.^) is given by 

p{X,S 0 {N,r))d\ = G''(iV,r)|A(A)|^5(r2-|Ap)dA. (1.11) 

Here, Cp(N, r) is a computable constant, S is the Dirac delta, and 

A(A) = n (1-12) 

The formula in (11.111) can be used—at least theoretically—to describe any statistic for these 
ensembles, such as spacing density, largest spacing, n-level correlation, n-level density, et. cetera. 
However, this is easier said than done. The analysis even to compute the marginals 

/ oo poo 

••• / p(A,5^(iV,r)) d\m+i---d\N (1.13) 

■OO J —oo 

was too cumbersome. The remainder of our results concern the empirical spectral density 

f(x,Sf(N,r)) = p'~»{(x),Sp{N,r)) (1.14) 

and its moments. 

For all spherical ensembles, we are able to express the moments in terms of the moments of 
the corresponding Gaussian ensembles, which allows us to prove a semicircular law in the limit. 
Explicitly, 

Theorem 1.4. Let (3 G {1,2,4}, and let q,r > 0. The moment of a spherical ensemble is 
related to the moment of a Gaussian ensemble by 

k 

mk{S 0 (N,r)) = ^ fi{N, q)), (1.15) 

where n = dimjj G^(iV, g) = /9('^) + N. Specializing to r = s/N, the eigenvalue density of 
Si 3 {N, s/N) has mean zero and variance one. The density functions f{x, S'^(iV, s/N)) converge in 
measure to fwig{x). 
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In the Hermitian case, we are able to give explicit formulas for the characteristic function and 
density function. 

Theorem 1.5. Let Ja be a Bessel function of the first kind (see chapter 9 of HASH ). The character¬ 
istic function equals 




r(^) 1 


N \rt 


j=0 


N 


(1.16) 


and the spectral density is 
f{x,SfiN,r)) = 


1 ^ f N 


E 


r(^ 


i) 


-+i 


(1.17) 


^iv2-2jY 2i0Fr i^^+j) ]\ \dxj 

- -/«[-) (1 - (; 

for a particular polynomial pn- 

The paper is organized as follows. We derive some necessary lemmas in ^21 We prove Theorem 
ll.3l in ^ Theorem 1 1.41 in ® and Theorem 1 1.5 1 in ^ We conclude by discussing spacing statistics 
in §0 


2. Preliminary Lemmas 

We calculate the joint and empirical spectral distributions of the spherical ensembles by relating 
them to the Gaussian ensembles, whose spectral distributions are known explicitly. 

Lemma 2.1. The ma^ of measure spaces A i—)■ defines a measure-preserving map from 

Gp{N,q) to Si 3 {N,r). In other words, if f E L^{Si 3 {N,r)), then 

[ f{.^)dsp{N,r)^ = [ f (dGp{N,q)A. (2.1) 

JSp{N,r) JGi 3 {N,q) Vl"^l / 

Proof Because the density function of a Gaussian ensemble depends only on the Frobenius norm, 
the mass of each ray from the origin is the same. Collapsing each such ray to the point on the 
Frobenius sphere it intersects gives the uniform distribution on the sphere. □ 

The next lemma will be needed in the proof of Theorem 1 1.41 

Lemma 2.2. Let G be a degree k homogeneous polynomial in n variables with complex coeffi¬ 
cients, B a positive definite symmetric matrix, \ ■ \ the standard 2-norm on M”, s/dffiBx the norm 
induced by B, and a > 0. Then 

L -(!)*:#) I ,2.2) 


9 

While this map is not dehned at the zero matrix, this is immaterial for our investigations below since the zero 
matrix is a single point of mass 0. 
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Proof. We first do the special case when B = I. Denote by dO the volume element on S'” For 
any function / : M —M so that G{x)f{\x\) is integrable, 


G{x)f{\x\)dx = 


Jo 


Jo 

poo 


G{r9)f{r)r^-^drd9 

1 

G{9)f{ry+'^-^drd9 


fc+n—1 


f{r)r 


If /(r) = e for a > 0, we have the r-integral is 


dr 


Ign-l 


G{9)d9. 


(2.3) 


f{r)r 


k+n-ldr = 


g-ar2/4^7+A:+n-l^^ 


The cases 7 = — fc and 7 = 0 give 


X 


^0 

1 fi 

2 Va 


2\a V2 



(2.4) 


fc+n 

= 1(1) “ + " 
2 V a 


G(0)d0, 

G0)d9. 


(2.5) 




Dividing these two equations gives the desired formula for the B = I case. 

The general case follows immediately from (12.51) . As B is positive definite, we have B = A~’'A 
for some full rank matrix A. The linear substitution y = Ax yields the more general formula, as 

G{X) -c.yBx)lA^^ ^ f 


y Bx) 


k/2 


I Ax I 


R- 1^1^ 

4 

a 
4 


cr \ 2 


r(i^) 

r(i) 


r(i±H), 

3. Proof of Theorem 1 1.3 1 


G {A-y det Ay dy 

G{x)e-<^^^^)yx. ( 2 . 6 ) 


□ 


The (order-symmetric) joint eigenvalue density for the N x N GOE / GUE / GSE ensembles are 
given by the formula 

p(X,a^{N,q)) = 7(iV,g)e-'’''l"l“/-‘|A(A)|'’ 
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(3.1) 














(see Theorem 3.1 of HRezll l. Here A = (Ai,..., Aat) is the veetor of eigenvalues, Cp{N, q) is a 
eomputable eonstant, and 

A(A) = (3.2) 

i<j 

In other words, for a symmetrie integrable funetion f of N real variables. 


<Gp{N,q) 


f{\{A))dGp{N,q)A. — / f{X)pN{X,G^{N,q))dX, 


(3.3) 


with p(A, Gi 3 {N, q)) as above. As by Lemma [ITI the spherieal ensemble measure is the pushfor- 
ward of the Gaussian ensemble measure under the map A —)■ ^^A, we have 


[ f{K^))dsp(N,r)A - f f fx f-—A\\ dGp{N,q)A 
JSpiN,r) dGp(N,q) \ \\^\ J J 

= f f (rTi^i^)] (^Gp{N,q)A 
JGp{N,q) Vl^l / 

^ j^jG^Ap{\G4N,q))d\. (3.4) 

Changing to spherieal eoordinates, we obtain 

/* p poo 

/ f{X{A))ds,^^,r)A= / f{e)p{te,Gp{N,q))rt^-^dtde, (3.5) 

JSp{N,r) Ar 5^-1 Jo 


where dO denotes the volume form of rS^ ^ (ealled dO earlier, in the r = 1 ease). As A 
is a homogeneous polynomial in the eigenvalues of degree (^), we have p{t9, G, q)) = 

G^(A, = G^(A, |A(6*)|^. Substituting into (13.51) . 

[ f{X{A))ds,iNM = [ r/(0)(G'(A,g)e-^'^^^*^/V(")|A(0)|^)rt^-'dfd0 

JSpiN,r) Jo ^ ^ 

= C'^(N,q)r(^P /(«) |A(0)|'’doj , 

(3.6) 


where (as usual) n = /9('^) + N. The substitution u = jdqrH"^/A immediately gives the value 
(I) for the first integral. Thus, 


JSp{N,r) 

where the normalizing eonstant 


f{X{A))ds,iN,r)A = G'^{N,r) 




f(e)\A{e)fde, 


')n—l 


0;t. i 

G';{N,r) = G'p{N,q)r -^T ' 


(3.7) 


(3.8) 


{jdqr'^Y'^'^ 

is independent of q. This equality (13.71) is equivalent to Theorem 1 1.3 1 The Dirae delta in Theorem 
ll.3l speeifies that the mass is eoneentrated on rS^~^. □ 
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4. Proof of Theorem 1 1.41 


Proof of Theorem lL^ We start by proving the moment expansion, equation (11.151) . Applied to the 
moments of the empirieal speetral density, Lemma [XT] implies 


mk{Sf}{N,r)) = 


1 

N 

1 

N 




TTiA’^)ds,iN,r)A 


Tr 


'Gp{N,q) 




dGfj{N,q)A 


J^Cp{N,qy 


lGp{N,q) 


|A|^ 


By Lemma IX 2 I we may rewrite the above as 


nik{Sp{N,r)) = ( ^ 


r(f) 1 
4 y r (^) N 


Cq{N,qy 


<Gp{N,q) 


Tr (A^) 


(4.1) 


(4.2) 


where n = dinijR Gp{N) = /9(^) + N (and we are using the faet that the norm of the Frobenius 
inner produet is a real inner product). Thus, 


mk{Si3{N,r)) 




mk{G(i{N,q)), 


(4.3) 


as desired. 

We now turn to the second part of the theorem, the convergence to the semi-circle. The eigen¬ 
value density of G^(A^, A^) has mean zero and variance 1-|- — A = l + o(l). The correct 

scaling for the spherical ensemble is to take r = \/N, as this will lead to eigenvalues of size on the 
order of 1. By (11.151) we find 


where n 


fB 2 \ 2 r f-) 

mt{Sf{N,VN)) = (^iVx/V j ^-^m^(Gj(JV.JV)) 

= (4.4) 

'^(^) + -^ = . In particular, Sj 3 {N, y/N) has mean zero and variance 

_ fyr2\ ( 1 A //iiV+(2-/3) 


/ /jiV+(2-/j) \ 
\ PN ) 


4 

(BN"^ 






iV(/3iV+(2-/3)) 



1 . 


8 


(4.5) 


























As —)■ oo, the asymptotic relation r(z + c)/V{z) ~ for \z\ —)■ oo implies the asymptotic 

equality 




Hn^- 

4 


k o \ ^ 

2 /yN^\ 2 


mk{Gfs{N,N)) 


mk{G0{N,N)). 


(4.6) 


In other words, if Gi denotes the Catalan number, then 


lim mfc(S'^(iV, \/iV)) = lim mk{Gfi{N,N)) = | 

" - " - '0 if /c IS odd. 


N^OO 


N^OO 


(4.7) 


which by the Method of Moments (see jBll |^) implies that, in the large N limit, the empirical 
spectral distribution for a spherical ensemble converges in measure to the semicircle. (See HRezH . 
Remark 2.2 for the fact that the limiting moments of Gaussian ensembles—and, more generally, 
Wigner ensembles—are Catalan numbers.) □ 


5. Proof of Theorem 1 1.5 1 

We first prove the expansion for the characteristic function in Theorem 1 1.5 1 and then we derive the 
claimed formula for the empirical spectral density. 

The moments of the GUE are well-known. The moment m 2 {Sp{N, r)) vanishes if k is odd. For 
k = 2i we have the following expression from Rarer and Zagier flHZH : 

m2iiG2iN , q )) = lg-^(2£-l)!!c(£,iV), (5.1) 

where c(£, N ) are power series coefficients for a particular rational function. 


and (2i — 1)!! = {2£ — l){2i — 3) ■ ■ • 3 ■ 1 is the 2(}^ moment of the standard normal. Substituting 
(15.11) into (11.151) yields 


m2t{S2{N,r)) = — ( — 


„2\ ^ 


r(i) 

r(| + (-) 


(2i-iyMi,N), 


and we may therefore write the characteristic function as 


<P{t,S2{N,r)) = ( ^ ( y 


£=0 


2\ ^ 


r(?) 

r(? + 4 




{ it ) 


21 


m 


r(t) 


y c(f, iV) ( 

7V'^r(f+ ^ ^ 4 


(5.3) 


(5.4) 
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For convenience, fix N and r, and set 


V'iv(w) 


1 

N 


OO 


E 


r(t) 


r (f + £) i\ 


c{i,N)u^, 


(5.5) 


so that '0AT (—S' 2 (iV, r)). The coefficients of this power series are a product of more 
familiar power series coefficients: 


CX) 


E 


r(i) 


r(f+ <)«' 



where is a confluent hypergeometric limit function, and by definition 


(5.6) 


= B{z) (5.7) 

i=0 

where B{z) is the rational function on the right hand side of (15.21) . From contour integration, we 
find 

We evaluate the contour integral by expanding both functions at 2 ; = 1, which is the location of the 
pole of B{z~^). The expansion of B is easy: 




N 


1 fN 


1 * = 


e=i 


z-1 


dz. 


The expansion of qFi is more involved but elementary: 
^ fn \ r (f) ^ „ 

qFi{—,Uz] = - iuzY 


e=o 

OO 


y ((z -1) +1)' 


e=o 

00 


r 1 °° 

OO -p /nX / OO 


j=0 \ £=j 

OO T~( / Ti \ f OO 


r if) I ^ r(;+j) 

OO p / n\ 

y —7—+ u^{z — iy. 


- ly 


(5.9) 


(5.10) 
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By the residue theorem, the only surviving terms after the integration are those with j 
Precisely, the integral in (15.81) evaluates to 


= -1. 


IpNiu) = 


2N 


N~l 

E 


j=o 

N-l 


j + iy r(2 + j)j! 






(5.11) 


and thus 


(p{t,S2iN,r)) = ipN[--irtf 


N 


N-l 

- V . 


r(t) ^ fn 1, ,2 

77+^>7(^^) 


1 ^ ^ 


(5.12) 


The expression in terms of Bessel functions in Theorem ll .5l follows immediately, using the formula 
(from HASH , p. 362) 


Ja{t) = 


my 


-qFi ( a + 1, 


(5.13) 


r(a + l)' 

This completes the proof of the first part of Theorem 1 1.51 

To compute the density function, we continue to use confluent hypergeometric limit functions 
because they have simple Fourier transforms. 

Lemma 5.1. The Fourier transform of qFi (a, is given by 

1 r ... _ / l .n\ . 1 /a-l\ 2^^ 

(1 - x^) 


271 


dt = 2 


„2^a-i 


(5.14) 


where the binomial coefficient is extended so that := r(^+i)r(t)-^+i) ’ meromorphic function 
of two complex variables w and z. 

Proof Fix a, and let w{t) = qFi (a, From the power series expansion of w, it’s easy to 

check that w satisfies the differential equation 

tw"{t) + (2q; — l)w'{t) + tw{t) = 0. (5.15) 

Taking the Fourier transform of both sides, and using the rules 

(5.16) 

= {ix)g{x), (5.17) 

we obtain the following differential equation for w: 

(1 — x^)w\x) + (2q; — 3)xw{x) = 0. (5.18) 

In other words, ^ (logt()(a;)) = {a — |)yE^, so w{x) = C{1 — x^)““t for some constant C. 

The formula 

1 /■°° ^ / 1 , lfa-1 

qFi ( a, ] dt = - 


C = — 

271 


follows from the formulas for the integrals of Bessel functions in Chapter 11 of HASH . 
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(5.19) 

□ 

















Figure 1 . Plots of the density f{x, 5'2(iV, y/N)) for = 4,6 and 8. 


It follows from Lemma [5TT] that 


1 

and 


n 


dt = 


1 ^ ^ 
i.2 


2i+i 


(±\ 


271 


n 




1 ^ ^ 
' \2 


2j 


dx J 


(l-x=) 




(5.20) 


2i+ir 

1 


+ J - 1 


dixjr 
f+j-l\ 




2i+lj..n-2 

From Fourier inversion we finally obtain 

/•OO 


dx j ^ 


2 2^ ^+i 

r — X W 


(5.21) 


/(x,^ 2 (iV,r)) = 


27r 


e-**^0(f,^2(A^,r))df 


— OO 

N-1 


N 


r(f) 


ly, , 

V2 + iyr(f+ j)j!2^+V 

r(f) 


n—2 


t+2-l 


^ ( 


n—3 I • 

2 2 \ +3 

r — X ^ 


iV-l 


iV 


2JY ^ \^j _|_ xy 2'?\/7rr + j) j! Vdxy/ 


d \ 


2i 


£) (5.22) 


and substituting n = N"^ eompletes the proof of Theorem 1 1.51 □ 

In Figures [U and m we plot the density /(x, S 2 {N, \/N)) for various values of N. (Note the 
polynomial expansions are valid for all N, though if is a perfeet square, then r = y/N is an 
integer.) The eonvergenee to the semieirele is quite rapid, and is similar to the eonvergenee to the 
semieirele found in HKKMSXII for bloek eireulant ensembles. 


6. Spacing Statistics 

We eonelude with a statistieal observation. The spaeing density of the SUE appears to approaeh 
the same large N limit as that of the GUE. This is not surprising, beeause spaeing densities are 
eonsidered to be very robust. Eor example, the thin ensemble of real symmetrie Toeplitz matriees 
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Figure 2. Plots of the density f{x, S 2 {N, y/N)) for = 9,16 and 25. 



Figure 3. Histograms showing spaeing statisties for the SUE, GUE, and the SUE 
and GUE together, respeetively. Eor both ensembles, 2000 matriees of size 100 x 
100 were sampled, and the 20 eonseeutive spaeings of the middle 21 eigenvalues 
were eomputed. 


appears to have GOE spaeings, even though its limiting empirieal density is far from semieireular 
HHML Spaeing statisties for the SUE and GUE are eompared in Eigure[3l 
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